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The cosmic expansion history tests the dynamics of the global evolution of the universe and its 
energy density contents, while the cosmic growth history tests the evolution of the inhomogeneous 
part of the energy density. Precision comparison of the two histories can distinguish the nature of the 
physics responsible for the accelerating cosmic expansion: an additional smooth component - dark 
energy - or a modification of the gravitational field equations. With the aid of a new fitting formula 
for linear perturbation growth accurate to 0.05-0.2%, we separate out the growth dependence on the 
expansion history and introduce a new growth index parameter 7 that quantifies the gravitational 
modification. 



I. INTRODUCTION 

Acceleration of the cosmic expansion reveals funda- 
mentally new physics missing from our picture of the 
universe, yet key for the understanding of the recent and 
present history and the fate of the universe. Further- 
more, this new physics tells us that our standard models 
of gravitation and particle physics may be woefully in- 
complete. The acceleration may lead us to insights about 
new high energy physics and the nature of the quantum 
vacuum, or about gravitation beyond Einstein's general 
relativity. Perhaps most exciting would be a signal that 
both are involved, providing clues to a theory of quantum 
gravity. 

The first scenario includes physical components such as 
the cosmological constant, dynamical scalar field models, 
or effective potentials from string theory. The second sce- 
nario includes extensions of the Einstein-Hilbert action, 
e.g. to higher derivative theories, scalar-tensor theories, 
generalized functions of the Ricci scalar, theories of su- 
pergravity or quantum gravity, and infrared modifica- 
tions of gravity such as from hidden spacetime dimen- 
sions. We can say that searching for the nature of the 
accelerating expansion is seeking to answer one or the 
other question: "Does nothing weigh something?" or "Is 
nowhere somewhere?" 

To distinguish the many different theoretical possibil- 
ities requires accurate observations of the cosmic expan- 
sion history, but even this will leave some degeneracies 
between explanations. Models with different physical ori- 
gins but the same global expansion properties could not 
be separated. Fortunately, the overall smooth character- 
istics of the universe are not the only observables. The 
energy density contents have evolved from the hot, dense, 
smooth state of the early universe to a relatively cool, 
diffuse, and in the case of matter, clustered state. While 
the first two properties are purely due to the expansion 
of the zeroth order, homogeneous universe, being qual- 
itatively kinematical, the last property arises from the 
perturbed, inhomogeneous universe, being intrinsically 
dynamical^]. The growth of large scale structure in the 
universe provides an important companion test, and the 
cosmic expansion history and growth history together 
provide discernment of the nature of the new accelerating 



physics. 

In jn]we discuss the expansion history and the effective 
equation of state entering the acceleration. The growth of 
linear perturbations in the matter component in a gener- 
alized cosmological model is reviewed in flllll The growth 
equation is extended in illVI to allow other theories of 
gravitation besides general relativity, and formal solu- 
tions given. For practical use in constraining models by 
observational data we introduce a highly accurate fitting 
formula in ij^and apply it to a braneworld gravity model 
and models with coupling between the matter and dark 
energy density. We present the conclusions in WII 

II. EXPANSION HISTORY 

The expansion history of the universe is a key quantity 
in cosmology, appearing directly in the metric in the form 
a(t). Kinematically, this is all that is needed to define 
distances and volumes (together with the spatial curva- 
ture constant k, which we take to be zero, though this 
does not affect the form of the following discussion) . To 
evaluate the distances for a specific cosmology, dynam- 
ics or equations of motion from the gravity theory are 
required, together with information on the energy den- 
sity contents. The expansion history, together with the 
amount of clustering matter and any interactions of the 
matter with other components, is the central ingredient 
as well in the growth of matter perturbations. 

The Friedmann expansion equation in terms of the 
Hubble parameter H = a/a is 

ff 2 = (87rG/3)5>, (1) 

i 

where we sum over all components of the energy density. 
Since we are especially interested in the matter compo- 
nent, e.g. since we are positive it exists and since we will 
later examine its growth into large scale structure, it is 
convenient to separate it out from the sum. Then in 
terms of dimensionless energy density we can write 

H 2 /H 2 = Q m a~ 3 + J2n»e a f» [1+wAa ' )]da ' /a ' (2) 
= n m a- 3 + 5H 2 /Hl (3) 
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where the set i' does not include matter, Q, m + Qj» = 
1, and w(a) is the equation of state of each component. 

Without imposing any physical interpretation on SH 2 
as actually being due to an energy density component 
as opposed to a modification of the Friedmann expan- 
sion equation, we can define an effective "acceleration 
physics" or "dark energy" equation of state 0] (cf. 0) 
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writing fi m (a) = Cl m a 3 /(H/H ) 2 . 

Volumes and distances are built up out of the confor- 
mal distance 



77(a) 



(da'/a')(a'H) 



dz'/H. 



(6) 



Models with the same expansion history will have the 
same distances and volumes. Note that formally we could 
obtain the same expansion history for two models by 
keeping their fl m the same and matching their w(a), or 
by allowing different f2 m and compensating for this in 
w(a). Since the latter case corresponds to misestimating 
the matter density rather than any new physics, we do 
not consider it further. 

While the definition of an effective equation of state in 
terms of the expansion history is powerful, allowing dif- 
ferent models to be talked about with a common language 
and treated in a model independent parameter space, this 
feature is also a bug. Measurements of the expansion 
history, through distances and volumes to arbitrary pre- 
cision, will not be able to distinguish different physical 
origins for the same expansion behavior. This is where 
the growth history comes in. 



III. GROWTH OF MATTER DENSITY 
PERTURBATIONS 



The universe has not remained homogeneous on all 
scales since its early, essentially smooth state. While the 
largest volumes can still be treated as homogeneous and 
isotropic Robertson- Walker universes, smaller scale evo- 
lution must take into account perturbations to the metric 
in the form of gravitational potentials. 

Note that recent speculation [j, |5j about the interac- 
tion of these potentials to affect significantly even the 
global expansion seems misplaced; investigation of a re- 
alistic inhomogeneous universe metric by Jacobs, Linder, 
& Wagoner ((J Q derived a Green function solution for 
the potential. This "post-Newtonian" solution corrects 
the Newton-Poisson equation and shows that no infrared 
divergences exist in the potential, rather a suppression 
as the Hubble scale is approached. The Appendix sum- 
marizes the effects. 



Using the perturbed equations of motion for the grav- 
ity theory, one can derive the growth of density pertur- 
bations. Concentrating on perturbations in the matter 
density S = 5p m / p m , assuming all other components are 
smooth, within general relativity the time evolution is 



6 + 2HS - 4irGp m 5 = 0. 



(7) 



The physical interpretation is very simple: the pertur- 
bations grow according to a source term involving the 
amount of matter able to cluster and are restricted by 
a friction term, or Hubble drag, arising from the expan- 
sion of the universe. General discussion of the physics 
dependence on the expansion rate is in 

It is convenient to study the growth evolution in terms 
of the expansion scale a or characteristic (e-fold) scale 
In a, rather than time t. Since the pure matter uni- 
verse solution has 5 ~ a, it is also useful in studying 
the dark energy to divide out this behavior and switch 
to the growth variable g — S/a. Finally, since we will be 
interested in modifications of gravity, we hereafter nor- 
malize G by Newton's constant, i.e. where G appears in 
equations it stands for G/GNewton- 

Denoting derivatives with respect to In a as primes, we 
have: 



4+i(lntf 2 )' 



9 



3+i(lnif 2 )'-^GO m (a) 
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-{l-w{a)]GQ w (a)g = 0, 
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0, (8) 
(9) 

(10) 

.9 = 0. (11) 



All these forms are equivalent. They clearly show that 
the Hubble drag is increased, and hence growth is sup- 
pressed, for an accelerating universe, as the deceleration 
parameter q — —ad/d 2 or w become more negative. We 
emphasize that they also demonstrate that within gen- 
eral relativity the linear theory growth factor depends 
purely on the expansion history, e.g. H(a) or w(a) or 
f2 m (a) or fl w (a) = 1 — Q m (a). So a discrepancy between 
the growth observed and that predicted based on an ob- 
served expansion history tests the theoretical framework 
and can point up modifications to the theory of gravity. 
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IV. GENERALIZATION TO GRAVITATIONAL 
MODIFICATIONS 

To study other theories of gravity we can consider 
a change to the effective Newton's constant G enter- 
ing the above equations (remember the G in the equa- 
tions really means G/GNewton), or more generally some 
non-zero right hand side. First we examine this as a 
generic change, and later treat a specific example within 
braneworld gravity. 

The deviation in G from its Newtonian value caused by 
some time variation in G can be viewed as a subset of a 
non-zero right hand side, since we may write a left hand 
side source term XG{a) as the usual X GNcwton and add 
a term X[GNcwton — G(a)] to the right hand side. Indeed, 
any difference between two cosmological models that only 
changes the source term, and keeps it linear in g, can be 
viewed as transforming the solution of the homogeneous 
differential equation for model 1 into a solution for model 
2 of the effective inhomogeneous differential equation. 

Using a Green function method one can obtain a formal 
solution 

g(ai,a) = g(ai,a)+ / du[Q(u) - Q{u)]g(a l ,u) 

x u 5 H(u)g(u)g(a) f dv g^ '(«) v~ 5 'tf" 1 ■(#?) 

J u 

where Q is the source term 1 . The barred quantities rep- 
resent some model 1 for which the solution is known (e.g. 
general relativity), and the integral gives the particular 
solution in the second model, for growth between any two 
scale factors (we can set a,; = to get the total growth 
up to some a). 

The solution can also be written as a recursion relation 

f a °° 
g{a,i,a) = g(a,i,a) + / dug(a z , u) ^ Kj(u, 0013) 

J °* i=l 

Ki + i(u,a) = / dx K\(u,x) Ki(x,a) (14) 

J u 

Kx(u,a) = [Q{u) - Q{u)]u 5 H(u)g(u)g(a) 

x j dvg^^v^H- 1 ^). (15) 

J u 

This is particularly useful when considering small per- 
turbations between models, e.g. when the gravitational 
coupling is slowly changing, as in the case of some scalar- 
tensor theories. (Retaining only the first order term, K\, 
is basically a Born approximation.) 

Another virtue of the Green function solution is the 
ability to see broad physical trends as the models change. 



1 Technically, Q is the source term divided by the growth vari- 
able g, and also multiplied by a~ 2 since Eq. 1121 uses a de- 
pendent variable of a rather than In a. For example Q = 
[(3/2) - (1/2) (In fil m (o))' - (3/2)Gfi m (a)]a" 2 in Eq. JTTJ. 



This follows the approach of 0, who considered the 
relation between distances as the cosmological model 
changed, including an analogous change in the theoreti- 
cal framework (there in terms of allowing a clumpy uni- 
verse). Here we consider the relation between growth 
factors. If Q > Q then Ki > and so g(di, a) > g(ai, a). 
Thus we have a criterion for when the growth will be 
stronger, or when it will be more suppressed. With the 
expansion history fixed, the criterion Q > Q simply be- 
comes G < G; i.e. if the effective gravitational coupling 
is stronger than Newton's constant then the growth is 
enhanced. For more elaborate modifications of gravity, 
a non-zero right hand side to the growth equation can 
contribute to Q as well, but the prescription above still 
applies. 



V. A NEW FITTING FORMULA FOR GROWTH 

The general growth solutions of the previous section 
are formal, and while we saw that they can present 
generic physical insights they are somewhat cumbersome 
for application to cosmological models. One might draw 
an analogy to trying to map the expansion history. While 
one can calculate the expansion history in a specific 
model, say from a high energy physics scalar field po- 
tential, this is inefficient for comparison of the observa- 
tions to many possible models. Instead a useful approach 
is a model independent one, using a parametrization of 
the expansion history, for example in terms of the equa- 
tion of state w{z) value and variation: w = w(z = 0) 
and w a = (—dw/da)\ z =i (this is also similar to the infla- 
tionary power spectrum index and tilt parameters). In 
this section we derive an analogous model independent 
parametrization of the growth history, putting it on equal 
footing with the cosmic expansion history. 

Rather than attempting to fit observations of growth 
history with an effective equation of state u> gr0 w (z) , it is 
better to render the physics appropriately: the expan- 
sion effects on the growth are described in terms of the 
standard expansion w(z), and the gravitational modifica- 
tions giving deviations from the expected growth history 
are treated as additional inputs. Again, in the standard 
framework the expansion history completely determines 
the growth history. Thus, we would like to write the 
growth history 17(a) as a function of an expansion history 
quantity plus a new, framework testing characteristic. 

Since the growth concerns matter density perturba- 
tions we take the expansion history in terms of the mat- 
ter density history f2 m (a). Two models with the same 
^jn(a) = Qm<i~ 3 / [H / Ho] 2 for all redshifts will have the 
same expansion history. So we look for a functional ex- 
pression g(Q m (a)). One that works superbly well, as 
both a highly accurate approximation to the exact solu- 
tion and as a simple characterization stimulating physical 
intuition, is 

/ n / dlna [n m .(a)' y — 1] /-in\ 

g(a) = e J o . (16) 
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Here 7 is our new parameter for the growth history, called 
the growth index, encompassing deviations in the theo- 
retical framework. Models with identical expansion his- 
tories but different gravitational theories will possess dif- 
ferent 7 parameters. 

A. Accuracy Tests 

First we establish the accuracy of the fitting formula, 
Eq. (|16JI . over a wide range of dark energy cosmologies. 
Note that 01 (also see HUE!) has found that a similar 
formula provides estimations of the normalized growth 
factor at the accuracy level of about 1%. However that 
approach normalized to the growth factor today (so it 
could not predict its value, and [2( showed that g(z)/g(0) 
varies by only a few percent innately between models) 
and fixed the grow th index. Furthering the pioneering 
work of [1(1 llll [l2j we can remove both those restrictions 
and obtain an order of magnitude better accuracy. 

In terms of the expansion history dark energy equation 
of state, within general relativity, we find excellent fits, 
to better than 0.2%, using 

7 = 0.55 + 0.05[l + w(z = 1)]. (17) 

Employing the value of the equation of state evaluated 
at z = 1 allows simple treatment of dynamical models 
where the equation of state varies with redshift, as it 
generically does. 

For the cosmological constant case, the fitting formulas 
Eqs. ijUfl) and l(T7l) reproduce the exact growth history for 
any redshift (including the total growth to the present) 
to better than 0.05% over the range £1 TO £ [0.22,1]! It 
remains accurate to better than 1% all the way down to 

n m = 0.01. 

Models with equation of state w — —0.8 (—0.5) have 
accurately fit growth histories to within 0.2% (0.4%) for 
Vt m G [0.2, 1]. A dynamical model such as SUGRA with 
w = -0.82, w a = 0.58 is fit to within 0.25%. Models 
with equations of state w < — 1 are similarly well approx- 
imated. When w = -1.2 (-1.5), the fit is good to 0.3% 
(0.5%). If we are willing to slightly modify the simplest 
fit of Eq. (H3J to 

7 = 0.55 + 0.02[l + w(z = 1)] forw<-l, (18) 

for the phantom models w < — 1, then we achieve an 
astonishing 0.05% accurac y fo r these fits. (Note that 
the fitting function of pH. Il4j|. also containing a single 
integral, is accurate to only 5% for models with w = —0.8 
or —1.2 and fl m = 0.3.) 

While impressive in accuracy, the growth function fit- 
ting formula's primary purpose is not a fit as such (the 
exact solution requires only solving a second order differ- 
ential equation), but rather its usefulness in physical in- 
tuition and in parametrizing modifications of the Einstein 
growth equation beyond the expansion behavior (just as 
w(z) parametrizes modifications of the Fricdmann expan- 
sion equation). The fitting function provides us access to 



the acceleration physics that exists beyond what the ex- 
pansion history sees. 

B. Example: Braneworld Gravity 

The growth of matter perturbations in gravitational 
theories beyond general relativity is not well developed. 
Here we consider one theory that has been shown to 
be self-consistent [H d Hi E Ojj , the DGP [H G3 
braneworld theory of gravity. In this theory gravity has 
infrared modifications due to spacetime possessing a large 
extra dimension (making our view a 4D brane within a 
5D bulk), causing a weakening of gravity on large scales 
approaching the Hubble scale. 

The expansion history for this braneworld theory fol- 
lows from the modified Friedmann equation, 

H 2 - H/r c = (8tt/3)p, (19) 

where r c = 7J ~ 1 /(1 — fl m ) is the crossover distance, re- 
lated to the 5D Planck mass. Equivalently the expansion 
history has an effective equation of state 

w(a) = -[l + n m (a)]~ 1 , (20) 

as noted elegantly by The braneworld expansion 

history can be well approximated by a simple scalar field 
model with wq = —0.78, w a = 0.32. Indeed these two 
very different physical origins for the acceleration agree 
in distance measurements to within 0.5% (0.01 mag) out 
to z = 2. 

Information from the growth history, however, as 
stated before can break this degeneracy in the nature 
of the acceleration physics. Comparing the braneworld 
model with a scalar field model with an identical expan- 
sion history shows deviations in the present growth factor 
of 7%. Figure n illustrates how the growth history de- 
pends on both expansion history and the gravitational 
framework. Taking into account only the expansion his- 
tory in the growth equation, the braneworld and scalar 
field models appear to have the same growth history. 
However, proper treatment of the gravitational modifica- 
tions inherent in the braneworld scenario separate these 
models. This has been pointed out as well in [l7ll20l l2l| . 

In the linear power spectrum the deviation ranges from 
4% at z — 2 to 15% today. While a scalar field model that 
matched the modified growth of the braneworld model is 
possible, it in turn can be distinguished through the ex- 
pansion history. We see that expansion measurements 
and growth measurements work in important comple- 
mentarity to reveal the nature of the new physics. 

Figure|21demonstrates this synergy explicitly. Suppos- 
ing the universe was described by a braneworld model 
with Q m = 0.28, distance measurements of the quality 
of the proposed Supernova/ Acceleration Probe (SNAP; 
[l^]) supernova data set, together with a 0.7% measure- 
ment of the distance to the CMB last scattering sur- 
face from Planck [2^, would provide the constraints (at 
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FIG. 1: The growth history is shown for an extra dimen- 
sional braneworld model (long dashed, blue curve) and a 
quintessence model with wo — —0.78, w a = 0.32 (short 
dashed, red), having nearly identical expansion histories. 
When proper account is taken of the effects of altered gravity 
on the braneworld growth history (solid, black curve) this al- 
lows distinction of these models. The expansion history can 
in turn rule out a quintessence model degenerate with the 
solid curve. 



FIG. 2: Expansion history and growth history constraints on 
the dark energy equation of state parameters can test the the- 
oretical framework by looking for inconsistent results. The 
blue cross gives the best fit for the expansion history of a 
quintessence (Q) universe matching the braneworld (BW) sce- 
nario, but the red star gives the best fit for the growth history 
to a quintessence model, assuming general relativity. The 
black ellipse shows the constraints at 68% confidence level 
from next generation data composed of SNAP supernovae 
data and Planck CMB last scattering distance measurement. 



68% confidence level, marginalizing over other parame- 
ters such as Q m ) in the dark energy equation of state 
parameter space shown by the ellipse. The best fit for 
growth measurements would concur with the solution 
using braneworld gravity equations, showing the consis- 
tency of the data with this model. However, if the growth 
equation employed general relativity, the best fit would 
lie at the red star, wo = —0.56, w a — 0, clearly inconsis- 
tent. The comparison of expansion and growth histories 
reveal a breakdown of the theoretical framework, this dis- 
crepancy alerting us to a possible modification of gravity 
(or experimental systematic errors). 

Of course if the measurements were too coarse and 
imprecise, we would not necessarily have noticed a sta- 
tistically significant discrepancy. The braneworld model 
expansion is compatible with an expansion history of a 
constant w = —0.71 model, to 1% in distance out to 
z = 1.7. So the expansion history measurements find a 
"distance" in equation of state space of Awo = 0.15 be- 
tween the effective scalar field model from the expansion 
history and that from the growth history. Conversely, 
the expansion history of (wo,w a ) = (—0.56,0) can be 
fit by (—0.63,0.32), so the "distance" from the expan- 
sion fit to the braneworld model of (—0.78, 0.32) is again 
Awq = 0.15. This suggests that for a 3<r detection of 
framework inconsistency we should strive for experiments 



that provide an uncertainty of <t(wq) < 0.05. 

Likewise one can estimate from the different orienta- 
tions of the expansion history and growth history con- 
straints in the wo — w a plane that the precision of mea- 
surements on w a should be a(w a ) < 0.2. This compari- 
son of growth to expansion provides one of the only ways 
of putting an absolute scale on the measurement pre- 
cision that should be striven for in experiments to re- 
veal the nature of the accelerating physics - a significant 
breakthrough (see also H3|). This is somewhat damp- 
ened by the realization that this scale is particular to the 
braneworld scenario. Note that the relative precisions 
between wq and w a obey the relation 

a(w') = a(w a )/2 « 2a(w ), (21) 

found in 0, H3- though that analysis was within the 
scalar field context. This relation signifies that a precise 
measurement of wq is of limited use without concomi- 
tant constraint on w a , since a sufficiently different w a 
can spoof wq. That is, the uncertainty in seeing a dis- 
crepancy will be dominated by the largest error among 
the two equation of state parameters. 

To move beyond a mere alarm that there is an incon- 
sistency, we need to employ the growth parametrization 
of Eq. I|16|) to obtain a quantitative measure of the devia- 
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tion from the growth behavior predicted by the expansion 
history measurements. We find that the fitting formula 
works for the braneworld scenario including gravitational 
modifications, using a growth index 7 = 0.68 (note that 
the pioneering paper of |l7j indicated the equivalent of 
7 = 2/3). In fact, the growth history using the fitting 
function Eq. (|16[) and 7 = 0.68 matches the exact solu- 
tion to within 0.2% (for fi ro > 0.2). 

The approximation of a single growth parameter be- 
yond the expansion history effects on the growth can 
be validated by asking what values of 7 as a function 
of redshift reproduce the exact solution. For the case 
n, n = 0.28, the (now) function 7(2:) ranges between 0.665 
at z = to 0.683 at z = 1 to 0.687 at high redshift. 
(The constancy of 7 with redshift holds even better for 
quintessence models.) This, as well as the excellent fit to 
the growth function, justifies the use of a single parame- 
ter 7, the growth index. 

With this model independent parametrization in hand, 
we can obtain quantitative measures of the deviations be- 
tween models, even those that involve gravitational mod- 
ifications. To the cosmological parameters f2 m , ion, and 
w a , we add the growth index 7 and can plot the resulting 
parameter estimation uncertainties, marginalizing over 
subsets of parameters. Figure illustrates an example. 

The growth index 7 that would fit the braneworld 
growth history is clearly distinct from the values allowed 
by a scalar field model that matches the expansion his- 
tory. The "distance" in 7 is 0.11; to attain the value 
7 = 0.68 would require, by extrapolation of Eq. (jT7Jl . 
w = +1.6! A 3cr distinction of the framework breaking 
would need a measurement with precision 17(7) = 0.037 
(marginalized over the other parameters). This corre- 
sponds roughly to a 2% measurement of the growth his- 
tory. Indeed, this is in good agreement with the results of 
Fig. ^ which showed growth deviations between the two 
models with identical expansion histories at up to the 7% 
level, so the same 3cr criterion leads to ~ 2% precision. 
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FIG. 3: While Fig. |5] showed that expansion history and 
growth history constraints on the dark energy equation of 
state parameters could test the theoretical framework by look- 
ing for inconsistent results, here we see quantitative mea- 
sures of framework breaking by gravitational modification of 
the growth index 7. The red star gives the best fit for a 
quintessence (Q) model matching the expansion history of the 
braneworld (BW) scenario, but the blue cross gives the true 
result for the braneworld growth history. The black ellipse 
shows the constraints at 68% confidence level from next gen- 
eration data composed of SNAP supernovae data and Planck 
CMB last scattering distance measurement, marginalized over 
the equation of state parameters Wq, w a . 



tions: 



pi = -3H(1 + Wi)pi + Ti(a, p m , Pde)- 



(22) 



C. Coupling of Matter and Dark Energy 

While the primary purpose of the formalism here is to 
test the gravitational framework, in the case of a physical 
dark energy there can enter microphysical effects. These 
can include spatial perturbations to the dark energy or 
coupling to the matter component. We leave the first of 
these to future work, but note that growth probes involv- 
ing correlations of large scale structure with the CMB 
might play a role (e.g. j2(| 153, HU). Here we consider 
whether the fitting formula and growth index approach 
remain valid in the presence of coupling. Without a mi- 
crophysical theory these are necessarily toy models, and 
we only consider the effects on matter growth, neglecting 
early universe or fifth force constraints. 

Interaction between matter and a dark energy compo- 
nent is treated through a coupling of the evolution equa- 



We have considered the cases T m = —Tde = ap m , apde, 
and aa n H, the slink y in flation model 29], and the undu- 
lant universe model 30] . Note the undulant universe is a 
special case of the slinky model, without coupling, and is 
ruled out by having a very low growth factor (go = 0.03). 

All these models follow the growth index formalism if 
the coupling is not too strong. As the coupling increases 
(e.g. if the dimensionless coupling a' = a/p m (0) > 0.5 
in the T m — aH case), this will start to break down 
because the equation of state of matter begins to devi- 
ate significantly from zero (plus in some cases the high 
redshift universe is not matter dominated). As a simple 
example, consider the decaying matter scenario where 
r m = —ap m . This was treated in detail in |3l| . and 
the matter equation of state is a/(3H) [12 • Generi- 
cally, a coupling of the form T,; = api will change the 
equation of state of component i, defined by the effec- 
tive conservation equation pi = —3H(1 + w e s)pi, from 
Wi to w c ff = Wi — a/(3H) (providing a way for a w > — 1 
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component to look phantom, w < — 1). 

Further research into the effects of coupling on growth 
of perturbations is underway [33| (also see H^D- Lack of 
a consistent microphysical theory is the major obstacle. 
For example, is the "new" energy density in a component 
distributed uniformly or in the same spatial distribution 
as the component it came from? Issues of evolution from 
that point by clumping or free streaming make rigorous 
calculation complicated. 



VI. CONCLUSIONS 

To reveal the physical origin of the acceleration of the 
universe, both probes of the expansion history (such as 
the distance-redshift relation) and of the growth history 
(such as weak gravitational lensing measurements involv- 
ing the mass power spectrum) are required. While the 
two types of probes in synergy give enhanced constraints 
on the effective dark energy equation of state, in compar- 
ison they can test the theoretical framework of cosmology 
and general relativity. 

The growth history of mass in the universe follows 
the source and friction term behaviors governed by the 
expansion. Deviations from this reveal a breakdown of 
the framework such as from modification of gravity. By 
rendering the growth function in a physically appropri- 
ate manner, separating the expansion effects from frame- 
work extensions, we presented here a new, physically in- 
tuitive and highly accurate (0.05-0.2%) fitting function, 
Eq. , for the linear growth of perturbations in gener- 
alized cosmologies. This allows model independent quan- 
tification of gravitational modifications in terms of a new 
parameter, the growth index 7. 

This research suggests a new paradigm for understand- 
ing the nature of the acceleration physics: accurate mea- 
surement of expansion and growth separately, for exam- 
ple through Type la supernovae and weak gravitational 
lensing. A useful, model independent, quantitative pa- 
rameter set was shown to be the equation of state value 
wo and variation w a and the growth index 7. In the 
specific worked case of comparing an extra dimensional 
braneworld scenario with scalar field physics in general 
relativity the desired measurement precisions should be 
of order a(w ) < 0.05, <r{w a ) < 0.2, (7(7) < 0.04. These 
should be technically feasible and should be within the 
reach of next generation experiments such as the Joint 
Dark Energy Mission. 

The formalism presented here has further applications 
for future investigation, such as seeing the effect of per- 
turbations in a physical dark energy component, cou- 
plings between dark energy and matter, and scalar-tensor 
gravity. To reveal the nature of the new physics responsi- 
ble for the universe-shaking acceleration, we will require 
a comprehensive suite of cosmological probes. The sig- 
nificance of the discoveries is so great that every robust 
method is needed to strengthen the accuracy, and the 
confidence in our understanding. With clear measure- 



ments of the cosmic expansion history and the cosmic 
growth history together we can learn if nothing weighs 
something, if nowhere is somewhere, or even more unex- 
pected insights. 
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APPENDIX: INHOMOGENEITIES AND COSMIC 
EXPANSION 

We have taken the cosmic growth history to not "back 
react" on the cosmic expansion history. That is, the 
global homogeneous expansion independent of the growth 
of matter structure is a valid treatment. This is a topic 
of great interest and comment; here we simply present a 
brief summary of the dependence of the metric on gravi- 
tational potentials and the lack of large contributions by 
gravitational potentials (in particular no infrared diver- 
gence) to the cosmic expansion. 

The approach taken by 0, Q is a straightforward cal- 
culation to obtain the metric of a realistically inhomo- 
geneous universe. In particular, it did not rely on any 
averaging procedure, rather a harmonic decomposition of 
the perturbations. The second key aspect was no a pri- 
ori assumption on the size of matter density fluctuations; 
rather it used a post-Newtonian parametrization, essen- 
tially a weak field, slow motion expansion. This followed 
work of Futamase [3^, |3|| |33 and can be traced back to 
the mean field theory, or two length scale, approach of 
Chandrasekhar |38[ . 

For potentials parametrized by an amplitude e 2 <C 1, 
and characteristic length scale ft, the slow motion or, 
more physically, the shear condition e 2 /n <C 1 applies. 
Violation of this condition leads to ray crossing in light 
propagation (see 0, H3) and eventually relativistically 
moving matter structures, contrary to observations of 
our universe. Landau & Lifshitz |40j pointed out that 
the dominant first order effect on the cosmic expansion 
entered at what they called pseudotensor order: e 4 / k 2 . 
Thus the shear condition ensures that the expansion 
is insignificantly affected, and conversely a significant 
back-reaction of inhomogeneities on the expansion would 
generically lead to visible anisotropies. 

However, here we concentrate on the post-Newtonian 
gravitational potentials, and modification of the Newton- 
Poisson equation relating the potentials to the matter 
density distribution. The general solution obtained by 
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was where a prime denotes a derivative with respect to the 

conformal time rj. 



(j>(r),x) 



4ir [ v du 



^ J no a '( u ) These expressions show that there is no divergence of 

d 3 " 3 ( )5 ( ~*)G( " ~*\ (All the potential or its derivatives in the presence of inhomo- 

ya [u) p[ ,y) y( ,rj, ,y), { ■ ) geneities. In contrast, the post-Newtonian Green func- 
tion solution, while matching the Newton-Poisson equa- 

plus an initial condition term. The Green function is tion on small scales, shows an exponential suppression 

j2 of the potential as one approaches horizon scales. These 

g(u, 77, x, y) = [a(u)/a(r])] [4ttC(u, rj)} limits are treated in detail in 7], and the physical prob- 

x g-ly-^l 2 /[4C(u,r;)] (A. 2) lem is shown to be closely analogous to the displacement 

pr] probability distribution for isotropic random walks, and 

C(u,rf) — (1/3) / dv(a/a'), (A. 3) for diffusion in a uniform medium. 

J u 



[1] Tegmark has written particularly clearly about this view 
of observational tests in terms of zeroth order and per- 
turbed universe properties. See, e.g., M. Tegmark, Sci- 
ence 296, 1427 (2002) astro-ph/0207199 [20] 
E.V. Linder & A. Jenkins, MNRAS 346, 573 (2003) 
astro-ph/0305286 [21 
U. Alam, V. Sah ni, T.D. Saini, A. S tarobinsky, MNRAS [22 
344, 1057 (2003) astro-ph/0303009 

E.W. Kolb, S. Matarrese, A. Notari, A. R iotto, Phys. [23 
Rev. D 71, 023524 (2005) hep-ph/0409038 [24 
E.W. Kolb, S . Matarrese^ A. Notari, A. Riotto, 
hep-th/0503117 [25 
M.W. Jacobs, E.V. Linder, R.V. Wagoner, Phys. Rev. D [26 
45, R3292 (1992) [27 
M.W. Jacobs, E.V. Linder, R.V. Wagoner, Phys. Rev. D 
48, 4623 (1993) [28 
E.V. Linder, A&A 206, 190 (1988) 

P. Schneider & A. Weiss, Ap. J. 327, 526 (1988) [29 
L. Wang & P.J. Steinhardt, Ap. J. 508, 483 (1998) [30 
D.H. Weinberg, to appear in New Astr. Reviews 
L. Amendola, C. Quercellini, E. Giallongo, MNRAS 357, [31 
429 (2005) astro-ph/0404599 [32 
D.J. Heath, MNRAS 179, 351 (1977) [33 
P.J.E. Peebles, Large-Scale Structure in the Universe [34 
(Princeton U. Press, Princeton, 1980) 

G. Dvali, G. Gabadadze, M. Porrati, Phys. Lett. B 485, [35 
208 (2000) hep-th/0005016 [36 
C. Deffayet, G . Dvali, G. "Gabada dze, Phys. Rev. D 65, [37 
044023 (2002) astro-ph/0105068 [38 
A. Lue, R. Scoccimarro, GD. Starkman, Phys. Rev. D [39 
69, 124015 (2004) astro-ph/0401515 [40 
G. Dvali, A. Gruzinov, M. Zaldarriaga, Phys. Rev. D 68, 
024012 (2003) hep-ph/0212069] 

G. Gabadadze, in From Fields to Strings: Cir- 



cumnavigating Theoretical Physics, eds. M. Shifman, 

A. Vainshtein, J. Wheater (World Scientific, 2004) 
hep-th/0408118 

Y.-S. Song"; Ph ys. Rev. D 71, 024026 (2005) 

astro-ph/0407489 

L. Knox , Y.-S. Song, J.A. Ty son, |astro-ph/0503644| 

SNAP - |http: / /snap Jbl.gov| ; 

Aldering TG. et al. 2004, |astr o-ph/0405232 
Planck - http://sci.esa.int/planck 

R.R. Caldwell & E.V. Linder, Phys. Rev. Lett, accepted 
astro-ph/0505494 

R.R. Caldwell fc E.V. Linder, in draf t 
L. Pogosian et al., astro-ph/0506396 

F. Giovi, C. Baccigalupi, F. Perrotta, Phys. Rev. D 71, 
103009 (2005) |astro-ph/0411702| 

B. Gold, Phys ! ReT D 71, 063522 (2005) 
astro-ph/0411376 

G. Barenboim & J. Lykken, astro-ph/0504090 

G. Barenboim, O. Mena, C. Quigg, Phys. Rev. D 71, 

063533 (2005) |astro-ph/0412010| 

M.S. Turner, Phys. Rev. D 31, 1212 (1985) 

E.V. Linder, A&A 206, 175 (1988) 

L. Barnes, M. Francis, GF. Lewis, E.V. Linder, in draft 

L. Amendola & C. Quercellini, Phys. Rev. Lett. 92, 

181102 (2004) astro-ph/0403019 

T. Futamase, Phys. Rev. Lett. 61, 2175 (1988) 

T. Futamase, MNRAS 237, 187 (1989) 

T. Futamase & M. Sasaki, Phys. Rev. D 40, 2502 (1989) 

S. Chandrasekhar, Rev. Mod. Phys. 15, 1 (1943) 

E.V. Linder, MNRAS 243, 362 (1990) 

L.D. Landau & E.M. Lifshitz, Classical Theory of Fields 

(Pergamon, Oxford, 1975) 



